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ABSTRACT 

We consider the power of a relativistic jet accelerated by the magnetic field of an accretion 
disc. It is found that the power extracted from the disc is mainly determined by the field 
strength and configuration of the field far from the disc. Comparing with the power extracted 
from a rotating black hole, we find that the jet power extracted from a disc can dominate over 
that from the rotating black hole. But in some cases, the jet power extracted from a rapidly 
rotating hole can be more important than that from the disc even if the poloidal field threading 
the hole is not significantly larger than that threading the inner edge of the disc. The results 
imply that the radio-loudness of quasars may be governed by its accretion rate which might be 
regulated by the central black hole mass. It is proposed that the different disc field generation 
mechanisms might be tested against observations of radio-loud quasars if their black hole 
masses are available. 
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1 INTRODUCTION 

The current most favoured models of powering active galactic 
nuclei (AGNs) involve gas accretion onto a massive black hole, 
though the details are still unclear. Relativistic jets have been ob- 
served in many radio-loud AGNs and are believed to be formed 
very close to the black holes. The spin energy of a black hole 
might be extracted to fuel the jet by magnetic fields supported by a 
surrounding accretion disc (Blandford & Znajek 1977; Blandford- 
Znajek process). This process has widely been believed to be 
the major mechanism that powers radio jets in AGNs (Begelman, 
Blandford & Rees 1984; Rees et al. 1982; Wilson & Colbert 1995; 
Moderski & Sikora 1996). However, Ghosh & Abramowicz (1997) 
doubted the importance of the Blandford-Znajek process. For a 
black hole of a given mass and angular momentum, the strength 
of the Blandford-Znajek process depends crucially on the strength 
of the poloidal field threading the horizon of the hole. They ar- 
gued that the strength of the field threading a black hole has been 
overestimated. The magnetic field threading a hole should be main- 
tained by the currents situated in the inner region of a surrounding 
accretion disc. Livio, Ogilvie & Pringle (1999; hereafter L99) re- 
investigated the problem and pointed out that even the calculations 
of Ghosh & Abramowicz (1997) have overestimated the power 
of the Blandford-Znajek process. Ghosh & Abramowicz (1997) 
have overestimated the strength of the large-scale field threading 
the inner region of an accretion disc, and then the power of the 
Blandford-Znajek process. They argued that the jet power extracted 
from an accretion disc dominates over the power extracted by the 
Blandford-Znajek process (L99). 

L99 estimated the maximal electromagnetic output from an 
accretion disc. In this case, the toroidal field component is of the 



same order of the poloidal field component at the disc surface due 
to the instability of a predominantly toroidal field (Biskamp 1993; 
Livio & Pringle 1992). Apart from the strength of the field thread- 
ing a disc, the acceleration of the jet is also governed by the mag- 
netic field configuration and the structure of the disc (Shu 1991; 
Ogilivie & Livio 1998, 2000; Cao & Spruit 2001; for relativistic 
jets near rotating black holes, see Koide, Shibata & Kudoh 1999; 
Koide et al. 2000). It is unclear how much power can actually be 
extracted from a magnetized accretion disc without considering its 
field configuration. In this work, we extend L99's work to explore 
some factors that affect the power of a relativistic jet extracted from 
an accretion disc. 



2 THE POWER OF A JET EXTRACTED FROM AN 
ACCRETION DISC 

For a relativistic jet accelerated by the magnetic field of the disc, 
the Alfven velocity is (Michel 1969; Camenzind 1986) 

Bp 

where and pa are the poloidal field strength and density of the 
jet at Alfven point, and jj is the Lorentz factor of the bulk mo- 
tion of the jet. Equation (1) is a good approximation, since the 
Alfven point is usually not very close to the rotating black hole 
(Camenzind 1986). The value of the Alfven velocity is of the order 
«a ~ RA^l(Rd), where Rd is the radius of the field footpoint at 
the disc surface, and Ra is the Alfven radius of the jet along the 
field line. 
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Mass and magnetic flux conservation along a magnetic field 
line requires 



mjet ^ Pava 



B A 



(2) 



where mj ct is the mass loss rate in the jet from unit surface area of 
the disc, and B p d is the poloidal field strength at the disc surface. 

Since the Alfven radius Ra depends on the global field con- 
figuration, not only the strength of the field at the disc surface, we 
assume the magnetic field far from the disc surface to be roughly 
self-similar: 



(3) 



where the self-similar index a > 1 (Blandford & Payne 1982; 
Lubow, Papaloizou & Pringle 1994). A self-similar index a = 4 
is adopted in the calculations of Lubow et al. (1994). 
The Lorentz factor of the jet is 



7j 



1-1^ 



(4) 



Combining Eqs (l)-(4), the mass loss rate in the jet from unit 
surface area of the disc is 



^pd / Rd^{Rd, 



4-7TC 



7f 



(7?-!) — 

The power of the jet accelerated by the magnetized disc is 



(5) 



L d = 27r / (7j - Vjrhjetc RddRd 



2nc 



Bpd 
4% 



Rd^l(Rd) 



/(«> jj)RddRd, 



where the function f(a, tj) is given by 

(7j - 1)7" 



/(«, 7j) 



(Tf 



For a given 7j , Eq. (6) reduces to 



L d = 27rc/(a, 7j ) / -JE2 



B p d ( Rd^(Rd) 



RddRd, 



(6) 



(7) 



(8) 



in the case of a = 1. The power of the jet Ld decreases with 
the self-similar index a, since the term RdQ(Rd)/c in Eq. (6) is 
always less than unity. Thus, Eq. (8) gives an upper limit on the 
power of a relativistic jet for given field strength at the disc surface 
and Lorentz factor 7j of the jet (tj S> 1; see Fig. 1). Note that Eq. 
(8) is almost same as Eq. (6) in L99 for the maximal electromag- 
netic output from the disc, except an additional term of f(a, 7 j) 
appears here. Here we have obtained a similar conclusion as that in 
L99, but in a rather different way. 

The toroidal field strength at the disc surface is 

47T m jct c 2 
Substituting Eq. (5) into Eq. (9), we have 



(9) 



B^d^f RdttjRd) 
Bpd V c 



f{a, 7j) 



(10) 



It is obvious that the ratio B,f,d /Bpd at the given radius Rd is deter- 
mined by the self-similar index a and the Lorentz factor jj, i.e., the 
ratio only depends on the jet properties, not on the disc properties. 



2.1 The magnetic field strength at the disc surface 

In order to complete the estimate on the power of the jet extracted 
from an accretion disc, we need to know the field strength at the 
disc surface. The origin of the ordered field that is assumed to 
thread the disc is still unclear. The strength of the field at the disc 
surface is usually assumed to scale with the pressure of the disc, 
as done in Ghosh & Abramowicz (1997). However, L99 have esti- 
mated the strength of the large-scale field threading the disc based 
on the results of some numerical simulations on dynamo mecha- 
nisms in accretion discs (Tout & Pringle 1996; Amitage 1998; Ro- 
manova et al. 1998). In this work, we mainly follow L99's estimate 
on the field strength at the disc surface. The physical implications 
of these two different estimates on the field strength at the disc sur- 
face will be discussed in Sects. 4 and 5. 

As L99, the strength of the magnetic field produced by dy- 
namo processes in the disc is given by 

^dvnamo 

^^~2H' (H) 
where W is the integrated shear stress of the disc, and H is the 
scale-height of the disc. For a relativistic accretion disc, the inte- 
grated shear stress is given by Eq. (5.6.14a) in Novikov & Thorne 
(1973, hereafter NT73). Equation (9) can be re-written as 



B, 



dynamo 



= 3.56 x 10 8 r" 3/4 m~ 1/2 J 4 _1 BL 1/2 gauss, 



where the dimensionless quantities are defined by 



R d „ GM M 

r d = r=r- , Rg = 

C 



(12) 



(13) 



Rg cr iviq 

and A, B and E are general relativistic correction factors defined 
inNT73. 

In standard accretion disc models, the angular velocity of the 
matter in the disc is usually very close to Keplerian velocity. For 
a relativistic accretion disc surrounding a rotating black hole, the 
Keplerian angular velocity is given by 

M 1 / 2 



0(r d ) 



3/2 



+ a 



(14) 



where a is dimensionless specific angular momentum of a rotating 
black hole. 

The large-scale field can be produced from the small-scale 
field created by dynamo processes as B(\) oc A -1 for the idealized 
case, where A is the length scale of the field (Tout & Pringle 1996; 
Romanova et al. 1998). The large-scale field threading the disc is 
related with the field produced by dynamo processes approximately 
by (L99) 



B 



P d 



f- 

\Rd 



* l 



B, 



dynamo ■ 



(15) 



The dimensionless scale-height of a disc h — H/Rd is in principle 
a function of Rd, and it reaches a maximal value in the inner region 
of the disc (Laor & Netzer 1989). We adopt the maximal value of 
H/Rd here. 

The scale-height of the disc H/Rd is given by (Laor & Netzer 
1989) 

H 



Rd 



26.345m acc Trf C2 , 



(16) 



where the coefficient C2 is defined in NT73. The dimensionless ac- 
cretion rate is defined by (Laor & Netzer 1989) 



Ah, 



1.626 x 10 15 m gs" 1 ' 



(17) 
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We use Eqs. (12) and (15), the power of the jet accelerated 
from a magnetized disc is available by integrating Eq. (6), if some 
parameters: m, rh acc , a, a and 7j are specified. 



3 THE RELATIVE IMPORTANCE OF THE 
BLANDFORD-ZNAJEK PROCESS 

As discussed in L99, the power extracted from a rotating black hole 
by the Blandford-Znajek process is determined by the hole mass 
M, the spin of the hole a, and the strength of the poloidal field 
threading the horizon of a rotating hole B p h' 



Lsz(max) 



r>2 

D ph n 2 2 

— — TvK h a c, 

47T 



(18) 



where is the radius of the black hole horizon. The strength of the 
poloidal field threading the hole is a crucial factor determining the 
power extracted from the hole. As argued in Ghosh & Abramowicz 
(1997) and L99, the field threading the black hole is generated by 
the currents situated in the disc rather than in the hole. Thus, it is 
realized that the field threading the hole should not be significantly 
larger than the field threading the inner edge of the disc Here we 
use a parameter £ to relate £? p h with B p ^(ri n ): 



B ph = (.Bpd(rin), 



(19) 



where f > 1, r- m is the radius of the inner edge of the disc. For a 
standard accretion disc, n n = r ms , where r ms is the radius of the 
minimal stable circular orbit. 

The Numerical simulations on the Blandford-Znajek process 
show that it can operate near its maximum power output (Komis- 
sarov 2001). Using Eqs. (6), (18) and (19), we can compare the jet 
power extracted from a magnetized disc with that by the Blandford- 
Znajek process. It is worth noting that the ratio of the jet power 
extracted from the disc to that from a rotating black hole is only 
determined by the parameters a, 7j and £. The relative importance 
of these two processes is independent of the mechanism for pro- 
ducing the disc field. 



4 RESULTS 

We plot the function f(a, jj ) varying with tj for different values of 
q in Fig. 1. The final velocity of a magnetically driven jet depends 
on the strength and global configuration of the field threading the 
disc, and the structure of the disc that provides boundary condi- 
tions for the jet (Shu 1991; Ogilivie & Livio 1998, 2000; Cao & 
Spruit 1994, 2001). In this work, we focus on relativistic jets which 
are widely observed in AGNs. To avoid being involved in the com- 
plexity of the jet acceleration problem, we use the Lorentz factor 
7j of a magnetically accelerated jet as a parameter to describe the 
problem. Now the problem is described by the parameters m, a, rh, 
a and 7j . 

The mass loss rate in the jet from unit surface area of the disc 
is plotted in Figs. 2 and 3 for the non-rotating and rapidly rotating 
black holes respectively. We plot the total mass loss rate in the jet 
and the power of the jet as functions of 7j for different values of a 
in Figs. 4 and 5 respectively. We calculate the ratio of the toroidal 
field strength to the poloidal one at the disc surface using Eq. (10), 
and the results are present in Fig. 6. 

Using Eq. (18), we can estimate the maximal power output 
from a rotating black hole by the Blandford-Znajek process. The 
power outputs as functions of a are present in Fig. 7 for different 
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Figure 1. The function f(a, 7j) varies with 7j for different values of a: 
a = l(solid line), 2(dotted), 3(dashed) and 4(dash-dotted). 



values of a and £. We plot the power outputs varying with accre- 
tion rate rh aC c based on L99's estimate of the disc field strength in 
Fig. 8. Finally, using Eqs. (6), (18) and (19), we plot the relative 
importance of these two power extract processes in the parameter 
space (a-a) in Fig. 9. 



4.1 The ratio of jet power to accretion luminosity 

The power of the jet extracted from a disc or a rotating black hole 
depends mainly on the field strength at the disc surface, since the 
field threading the hole is generated by the currents situated in the 
disc. The ratio of jet power to accretion luminosity is therefore re- 
lated with the mechanism for generating the disc field. 

In the case of the large-scale field being produced from the 
small-scale field created by dynamo processes, the large-scale field 
strength is estimated by Eq. (15). As we know, H/R^ oc m acc and 
L acc oc m acc m, then we obtain a linear relation between m acc and 
ijet/^aco by using Eqs. (6) (see Fig. 8): 



Ljet/Lacc OC 77Za.cc, 



(20) 



which is independent of the black hole mass. Note that the strength 
of the field threading the hole is related with that threading the disc 
(see Eqs. (18) and (19)), the relation (20) is also valid for jet power 
extracted from the rotating hole. 

The similar analysis can be performed for the case that the 
magnetic field strength scales with the pressure of the disc. The 
strength of the field at the disc surface is 



B 



pd (r d ) oc r d 



3/4 



m-^A^BE 1 ' 2 , 



(21) 



which is in dependent of the accretion rate (NT73). 

From Eqs. (6), (18) and (21), we obtain a relation (Ghosh & 
Abramowicz 1997) 



Bjet /La 



(22) 



which holds for a jet accelerated either from a disc or a rotating 
hole. It does not depend on the black hole mass either. 
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Figure 2. The mass loss rate in the jet from unit surface area of the disc 
surrounding a non-rotating black hole for different values of a. The hole 
mass m = 10 9 , j; = 10 and m aC c = 0.05 are adopted. 



Figure 4. The total mass loss rate in the jet as functions of 7j for different 
values of a for the field. The black hole spin parameter a = 0.95 and 
m acc = 0.05 are adopted. 
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Figure 3. Same as Fig. 2, but for a rapidly rotating black hole: a = 0.95. 
5 DISCUSSION OF THE RESULTS 

Figure 1 shows that the function f(a, jj) varies slowly with 7j ex- 
cept for 7 — > 1, The function f(a, 7j) is almost a constant while 
7j > 5. From Eq. (6), it is then found that the power of a relativistic 
jet accelerated by a magnetized accretion disc is insensitive to its 
Lorentz factor 7j . This is confirmed by the results plotted in Fig. 5. 
The jet power varies with jj in a range less than 10 % for 7j > 5. 
It implies that relativistic jets with different Lorentz factors 7j ac- 
celerated by magnetized discs have almost same jet power, if all 
other physical parameters are fixed. We can therefore adopt a typi- 
cal value of 7j in our calculations, which will not change the main 
results on jet power, though the terminal velocity of the magneti- 
cally driven jet could be given in principle only after the structure 
of the disc and its field strength and configuration are supplied. This 



Figure 5. The power of the jet accelerated by a magnetized disc as func- 
tions of 7j for different values of a. The left panel is for a non-rotating black 
hole, while the right panel is for a rapidly rotating black hole: a = 0.95. 
The accretion rate m aC c = 0.05 is adopted. 

has simplified our present investigation, since we are interested in 
the power of relativistic jets extracted from discs. 

The power of a jet decreases with the self-similar index a (see 
Fig. 5). A large a implies that the poloidal field strength decays 
rapidly with radius along a field line. From Fig. 4, we see that less 
mass is loaded into the jet in the case of a larger a, and the power 
extracted from the disc is then reduced. The power of the jet is 
sensitive to the field configuration far from the disc surface. 

We find that the mass loss rate in the jet from unit surface area 
of the disc decreases with radius (see Figs. 2 and 3). Almost all 
mass carried by the jet is from the region of the disc: R c \ < 5 Ri n , 
i.e., jets are mainly accelerated from the inner region of the disc 
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Figure 6. The ratio of the toroidal field strength to the poloidal one at the 
disc surface as functions of 7j for different values of a at ij d = 2R ln . The 
hole spin parameter a = 0.95 is adopted. 



Figure 8. The power outputs as functions of m acc for different values of 
a and £. The solid lines represent the power output from an accretion disc, 
while the dotted lines represent the power output of the Blandford-Znajek 
process. The black hole spin parameter a = 0.95 is adopted. 



2, — 4 = — 



i-3 



o 





Figure 7. The power outputs as functions of a for different values of a and 
f . The solid lines represent the power output from an accretion disc, while 
the dotted lines represent the power output of the Blandford-Znajek process. 
The accretion rate m acc = 0.05 is adopted. 



Figure 9. The relative importance of two power extract processes in param- 
eter space (a-a) for different values of £ = l(solid line), 2(dotted) and 
3(dashed). 7j = 10 is adopted. 



where most gravitational energy of accretion matter is released. The 
total mass loss rate in a relativistic jet can be neglected compared 
with the accretion rate of the disc (usually Mj c t < 10 _2 Af acc , see 
Fig. 4). A jet with high mass loss rate results in a low terminal 
velocity of the jet and wound-up field lines (Figs. 4 and 6), which 
is consistent with the results of Cao & Spruit (1994). As pointed 
out in Sect. 2, our calculation of the jet power extracted from a disc 
can reproduce L99's result while a = 1 is adopted (see Eq. (8)). In 
L99's calculation, B^,d = B p d at the disc surface is assumed, so it 
is not surprised to find that B^d — B p d in the case of a — 1 in our 
calculations (see Fig. 6). For a large a, the ratio B^d / B p d is low 
(compare the results in Fig. 6 for different values of a). 



In our present calculations, the field configuration far from the 
disc is described by the self-similar index a. In the case of a = 1, 
it is found in Fig. 7 that the power extracted from the disc always 
dominates over the maximal power of the Blandford-Znajek pro- 
cess, if £ < 3, i.e., the strength of the poloidal field threading the 
hole is less than three times of that threading the inner edge of the 
disc. If £ = 1 is assumed, i.e., B p h — B p d(Rin), Ld is an order 
of magnitude larger than Lbz even for an extreme Kerr black hole. 
Thus the conclusion given by L99 that the maximal power extracted 
from a disc dominates over that from a rapidly rotating black hole 
is confirmed, while our calculations are not limited to the extreme 
case a = 1. We present the results in Fig. 7 to compare the power 
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outputs by these two processes for different values of a and £. It 
is found that the Blandford-Znajek process would be important for 
a large a or/and From Eqs. (6) and (18), we find that this con- 
clusion is independent of the mechanism for producing disc fields. 
The conclusion is also independent of the hole mass m and accre- 
tion rate rh acc , since the hole field strength is only governed by the 
disc field strength. 

If the large-scale field can be produced from the small-scale 
field created by dynamo processes as done by L99, we find that 
linear relations are present: Ljet/iacc oc rhacc, for a jet accel- 
erated either from a disc or a rotating hole, (see Fig. 8 and Eq. 
(20)). This relation is quite different from the situation for the 
field that is assumed to scale with the pressure in the disc, where 
Lj e t I Lace oc m~c is present. 

The ratio Ljet/Lacc can be related with an observational 
quantity of quasars: radio-loudness R, which is defined as R = 
/„(5GHz)//„(4400A). The fact that the ratio Lj c t/Lacc is inde- 
pendent of the black hole mass may imply that the radio-loudness R 
is not related directly with the black hole mass. The radio-loudness 
R may probably governed by the accretion rate. The accretion rate 
is then regulated by the central black hole mass in some way (Yi 
1996; Salucci et al. 1999; Haiman & Menou 2000; Kauffmann & 
Haehnelt 2000). The relations found between radio-loudness and 
the central black hole masses in quasars (McLure et al. 1999; Laor 
2000; Gu, Cao & Jiang 2001; Ho 2002) may reflect intrinsic re- 
lations between radio-loudness and accretion rate which is regu- 
lated by the black hole mass (Boroson 2002). The different rela- 
tions between Lj e t / Lace and rh aC c are expected for different mech- 
anisms of producing the disc field. The different disc field gener- 
ation mechanisms can therefore be tested against observations of 
radio-loud quasars if their black hole masses are available. How- 
ever, whether the jet power is extracted from a disc or a rotating 
black hole cannot be tested in this way, since the calculations pre- 
dict the same relation between Lj ct / Lace and rh aC c for these two 
jet acceleration mechanisms. 

In Fig. 9, we compare the relative importance of these two 
power extract processes in the parameter space (a — a). Only 
for those cases the strength of the poloidal field of the disc de- 
cays slowly with radius along the field line (small a), the power 
extracted from the disc dominates over that extracted by the 
Blandford-Znajek process even for a rapidly rotating hole. In the 
right-upper corner of the figure, i.e., the case with a large a and a, 
the power extracted from the disc is greater than the power of the 
Blandford-Znajek process. Otherwise, the power of the Blandford- 
Znajek process can be neglected compared with that from the disc. 

The acceleration of the jet is governed by the structure of the 
disc and the strength and configuration of the field threading the 
disc. The physical properties of the jet, such as the power and ter- 
minal velocity of the jet, can be determined if all these physical 
factors of the disc and its field are known. Instead of solving a set 
of MHD equations describing the disc-jet system, a self-similar in- 
dex a is employed to describe the field far above the disc surface 
in present work. The physics of how the jet acceleration is gov- 
erned by the detailed structure of the disc near a rapidly rotating 
black hole has not been included in this work. It would be the main 
limitation of the present work. 
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